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This paper shows how to analyze plasmonic metal-insulator-metal waveguides using the fuU modal structure of 
these guides. The analysis applies to all frequencies, particularly including the near infrared and visible spectrum, 
and to a wide range of sizes, including nanometallic structures. We use the approach here specifically to analyze 
waveguide junctions. We show that the full modal structure of the metal-insulator-metal (mim) waveguides — 
which consists of real and complex discrete eigenvalue spectra, as well as the continuous spectrum — forms a 
complete basis set. We provide the derivation of these modes using the techniques developed for Sturm-Liouville 
and generalized eigenvalue equations. We demonstrate the need to include all parts of the spectrum to have a 
complete set of basis vectors to describe scattering within mim waveguides with the mode-matching technique. 
We numerically compare the mode-matching formulation with finite-diflierence frequency- domain analysis and 
find very good agreement between the two for modal scattering at symmetric mim waveguide junctions. We 
touch upon the similarities between the underlying mathematical structure of the mim waveguide and the 77 
symmetric quantum mechanical pseudo-Hermitian Hamiltonians. The rich set of modes that the mim waveguide 
supports forms a canonical example against which other more complicated geometries can be compared. Our 
work here encompasses the microwave results, but extends also to waveguides with real metals even at infrared 
and optical frequencies. 
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I. INTRODUCTION 

Waveguides have long been used to controllably direct en- 
ergy flow between different points in space. Understanding the 
way waves propagate in waveguides led to a multitude of cre- 
ative designs — all the way from the pipe organ to Ught switches 
used in fiber optic communications. In optics, recently there 
has been a growing interest in making use of the dielectric 
properties of metals to guide electromagnetic energy by using 
sub-wavelength sized designs that work in the infrared and the 
visible bands of the spectrum. One of the motivations for do- 
ing photonic research using metals is to find the means to inte- 
grate electronic devices with sizes of tens of nanometers with 
the relatively much larger optical components — so that some 
of the electrons used in the communication channels between 
electrical circuitry can be replaced by photons for faster and 
cooler operation'. 

Whereas the use of metals for directing electromagnetic en- 
ergy is relatively new in optics, sub-wavelength guiding of light 
by metals is the norm in the microwave domain. Even though 
the permittivity of metals can be large in magnitude at both 
microwave and optical frequencies, the characteristics of the 
permittivity are quite different in the two frequency regimes. 

In the microwave regime, electrons go through multiple col- 
lisions with the ions of the lattice during an electromagnetic 
cycle according to the phenomenological Drude model of elec- 
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trons. Therefore, the electron movement is a drift motion 
where the velocity of electrons is proportional to the applied 
field strength (Ref. |2j Chp. 1). As a result, the induced dipole 
moment density and hence the permittivity is a large, negativ^l 
imaginary number On the other hand, at optical frequencies 
below the plasma oscillation frequency, electrons go through a 
negligible number of collisions during an electromagnetic cy- 
cle and this time acceleration of electrons is proportional to the 
appUed field strength which then results in a permittivity that 
can be substantially a large real negative number. Above the 
plasma frequency, the induced dipole moment density is very 
low and the permittivity is predominantly a positive real num- 
ber less than one (Ref.j4j Chp. 9). 

The dielectric slab and the parallel plate (i.e. consisting of 
two parallel perfectly conducting metal plates) waveguides are 
the two canonical examples of waveguiding theory. If we have 
a layered metal-insulator-metal (mim) geometry, it is possible 
to smoothly transition from the dielectric slab to the parallel 
plate waveguide by reducing the frequency of operation, and 
therefore varying the metallic permittivity e^. 

At frequencies above the plasma frequency, the metal has 
a permittivity < 1 whereas the insulator has e, > 1. We 
illustrate the geometry in the inset of Fig. [1] 

The physical modes that the dielectric slab waveguide sup- 
ports fall into two sets: guided modes and radiation modes 
(Ref. |5] Chp. 1). Guided modes consist of a countable, finite 
set, i.e. there is only a finite number of discrete guided modes. 
Radiation modes consist of a non-countable, infinite set, i.e. 
they form a continuum. The combination of these two sets of 
modes form a complete and orthogonal basis set. 

Now suppose that we change our operation frequency to one 
which is very close to the dc limit where is an arbitrarily 
large, negative, imaginary number. In this limit, we can ap- 
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proximate the metal as a perfect electric conductor (pec) where 
|e„| -i- oo. Such an approximation then gives us the parallel 
plate waveguide of the microwave domain. Unlike the dielec- 
tric slab, the parallel plate geometry is bounded in the trans- 
verse dimension — fields are not allowed to penetrate into the 
PEC. There are infinitely many discrete modes of the parallel 
plate, all of which have sinusoidal shapes, and there are no con- 
tinuum modes. The collection of the infinitely many discrete 
modes forms a complete orthogonal basis set. 

In this paper we will investigate the modal structure of the 
two dimensional mim waveguide in the infrared regime where 
Cm is primarily a large, negative real number The geometry of 
the MIM waveguide is exactly the same as the one in the inset 
of Fig. |i] The only difference between the parallel plate, mim 
and the dielectric slab waveguides is in the numerical value of 
£„, which depends on the frequency of operation. 

There have been numerous s tudie s of the mi m waveguide 
in the liteiiitm^^^ismi^mimL^SM^^sm^^iisJm. The fact that 
light can be guided within a deep subwavelength volume over 
a very wide range of wavelengths is one of the primary reasons 
why the mim geometry has attracted so much attention. The 
full set of modes that the mim waveguide supports — real and 
complex discrete modes as well as a continuous set of modes — 
has only very recently been published^' . For other geometries, 
it has been shown that, in general, waveguides support real, 
complex and continuous sets of modejS F^P^M^sl jj^ ^jjjg work, 
we will provide the detailed mathematical framework to ana- 
lyze the modal structure of the mim waveguide and emphasize 
how it is a hybrid between the parallel plate and the dielectric 
slab waveguides. 

Operator theory will be the basis of the mathematical tool 
set with which we will analyze mim waveguides. In Section 
II, we will introduce the notation and make some definitions 
pertaining to the operators in infinite dimensional spaces. In 
Section iii we will derive the discrete and continuum modes 
supported by the mim waveguide and show that the underlying 
operators are pseudo-Hermitian. In Section iv we will demon- 
strate that the modes we report form a complete basis set via 
example calculations using the mode-matching technique. In 
Section v we will discuss our results and underline some of 
the relevant developments in mathematics from both quantum 
mechanics and microwave theory with the hope of expanding 
our tools of analysis. Lastly, in Section vi we will draw our 
conclusions. 



II. SOME DEFINITIONS 

Throughout the paper, we will be using nomenclature from 
operator theory. In this section, we will define the terminology 
and introduce the notation which will be used in the following 
sections^. The reader well versed in operator theory can di- 
rectly skip to the next section. 

A linear vector space is a space which is closed under the 
operations of addition and of multiplication by a scalar. We 
will call the elements of the space vectors. Spaces need not be 
finite dimensional — infinite dimensional vector spaces are also 
possible. For instance, the collection of all square integrable 



DC Microwave Infrared Visible Ultraviolet 
I 1 1 1 1 1 — I >co 



1/T CO 

p 

FIG. i: The change in the metallic permittivity (e^ ) as the frequency 
of operation (o)) is varied leads to an evolution from the parallel plate 
waveguide at low frequencies to the dielectric slab waveguide at high 
frequencies for the two dimensional metal-insulator-metal geometry. 
We are assuming that e,„ is a Drude model metal where t is the aver- 
age time between collisions among the ions of the lattice and the free 
electrons, copis the plasma oscUlation frequency. 



functions f{x) defined on an interval a < x < b forms an 
infinite dimensional vector space. 

The inner product is a scalar valued function of two vectors 
/ and g, written {f\g) with the following properties 

{f\g) = {g\fr 

{aj + a2g\h) = a*{f\h) + a*{g\h) 
(/|/)>0 if/+0. 

Here (•)* denotes complex conjugation, a{i,2} are arbitrary 
complex numbers and f,g,h denote arbitrary members of the 
linear vector space §. 

For the infinite dimensional vector space of square inte- 
grable functions one possible definition of the inner product 
is 

(f\9)= f r{x)g{x)Ax. (1) 

J a 

A hnear vector space with an inner product is called an inner 
product space. In such spaces the norm of a vector f is defined 
as 

ll/ll = vW)- (2) 

This is also known as the norm, to denote square integra- 
bility in the sense of Lebesgue. By using the norm of a vector 
space, we can define the distance between its vectors / and g 
as d{f,g) = 11/ - g\\ which is always nonzero if / + g. Here, 
d{f, g) is called the metric — the measure of distance between 
vectors — of the inner product space. Suppose that 3^ and 9 are 
two subsets of the inner product space S and that 3^ is also a 
subset of 9, i.e. 3^ c 9 c 8. 3^ is said to be dense in 9, if for 
each (J 6 9 and £ > 0, there exists an element / e 3^ where 
d{f,g) <£ (Ref |3i] pp. 94-95)- 

A vector space § is complete if all converging sequences of 
vectors /„ (x) converge to an element / e S. An inner product 
space which is complete when using the norm defined by 
Q is called a Hilbert space. 

An operator L is a mapping that assigns to a vector / in a 
linear vector space Si another vector in a different vector space 
§2 which we denote by L/ (most often Si = S2). An operator 
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is linear if \-{a\f + aig) = ociLf + a2Lg for arbitrary scalars 
a{i_2} and vectors /, g. The domain of an operator L is the set of 
vectors / for which the mapping L/ is defined. The range of an 
operator L is the set of vectors g = Lf for all possible values of / 
in the domain of L. A linear operator is bounded if its domain is 
the entire linear space S of vectors / and if there exists a single 
constant C such that ||L/|| < C||/||. Otherwise the operator 
is unbounded. The differential operator is a classical example 
of an unbounded operator (Ref |32| pp. 93-94). L is positive 
(negative) definite if (/|L/) > ((/|L/) < o) for all possible/. 
Otherwise, L is indefinite. 

A linear bounded operator is said to be the adjoint of L if 
for all / and g'm$ the condition {g\\-f) = {L^ g\f) is satisfied. 
If L = then L is said to be self-adjoint. If the operator L 
is unbounded, then the equality {g\\-f) = {V g\f) defines a 
formal self-adjoint (Ref 31 , Sec. 3.4.1). 

Suppose we have a set of orthonormal vectors {/„} which 
span the Hilbert space "K. Then, we can expand any vector g e 
g = T.„{fn\g)fn- Similarly, any linear bounded operator 
L acting on g results in 

L^= E(/"l^)'-/« = Y,ifn\9){fm\^fn)fm 

n n,m 

= T.fm{fm\\-f„){fn\g) 

where we expanded L/„ in terms of the basis {/„ } to get to the 
last line. Once we choose a complete orthonormal basis set, 
we can describe the action of L on any vector g by the product 
of an infinite dimensional matrix with elements (/m | L/„ ) and 
an infinite dimensional vector with elements (/„ \g) — a gener- 
alization of regular matrix multiplication. The infinite dimen- 
sional matrix is called the representation of L in {/„}. If the 
matrix for L is diagonal, then we call that the spectral represen- 
tation (Ref 33. p. 110). 

The spectral representation for an operator L depends on the 
study of the inverse of the operator L - A, which we will denote 
by (L- A)"\ for all complex values of A (Ref [33] p. 125). Let the 
domain and range of L be denoted by Di and 3?^. The point 
(discrete) spectrum is the set of A for which (L - A)~^ does not 
exist. The continuous spectrum is the collection of A for which 
(L - A)"^ exists and is defined on a set dense in IRi, but for 
which it is unbounded. The residual spectrum is the collection 
of A for which ( L - A)"^ exists (it may or may not be bounded), 
but for which it is not defined on a set dense in "Ri. The spec- 
trum of L consists of values of A which belong to either the 
point, continuous or the residual spectrum (Ref 32 , p. 371 and 
Ref |34j p. 21). We summarized the taxonomy of the spectrum 
in Fig. [2] 

III. SPECTRUM 

After having defined the necessary terminology, in this sec- 
tion we will derive the modal structure (spectrum) of the mim 
waveguide. We will specifically focus on the even modes of the 
waveguide, for which the transverse magnetic (tm) field com- 
ponent is an even function of the transverse coordinate, x. The 
reason why we focus on even modes is that we will be analyzing 
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FIG. 2: Classification of the spectrum for an operator in a general 
space. Point spectrum is sometimes called the discrete spectrum. 



the scattering of the main, even mode of the mim waveguide — 
which is also a tm mode — off of a symmetric junction with a 
different sized mim waveguide. Due to the symmetry of the 
problem at hand, even modes will be sufficient. We could also 
solve for the case of the odd modes by a similar approach, but 
we omit that explicit solution for reasons of space. Evenness 
of the function is achieved by putting a fictitious perfect elec- 
tric conductor (pec) at the x = plane of the waveguide, which 
forces the tangential electric field to be an odd function, 
and the magnetic field Hy to be an even function of x. In other 
words, the modes of this fictitious waveguide with the pec at 
X = are mathematically the same as the even modes of the 
actual waveguide of interest, and so we will work with this hy- 
pothetical waveguide for our mathematics. The geometry is as 
shown in Fig. [3] £,„ refers to the permittivity of the metal re- 
gion and e, of the insulator region. At infrared frequencies, Cm 
is a complex number with a large, negative real part and a rela- 
tively small imaginary part (the sign of which is determined by 
the time convention used, being negative for an exp(ia)f) time 
dependence). 

Let us begin with Maxwell's equations for fields that have an 
exp(ia)t) time dependence. 

V X E(r) = -ia)u(x)H(r) 

(^) 

VxH(r)= ia)e(x)E(r). 




r V 

FIG. 3: (Color online) Geometry for the even modes of the mim wave- 
guide. The X = plane contains a fictitious perfect electric conductor 
to simplify the problem when dealing only with even tm modes of the 
guide. This fictitious mim waveguide is equivalent to an actual guide 
with an insulator thickness of la. The inset shows the equivalent sym- 
metric junction of two mim waveguides. The dashed line in the inset 
is the plane of symmetry, which is where the fictitious pec layer is 
introduced. 
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The MiM waveguide is a two dimensional structure which does 
not have any variation in the y direction. Therefore, we can 
eUminate all the derivatives with respect to y in Maxwell's 
equations. Furthermore, our study will be based on the tm 
modes which only have the H^, Ex and field components. 
Also, the uniformity of the waveguide in the z direction leads 
to exp(-ifczz) as the space dependence in z by using the sepa- 
ration of variables technique for differential equations (fc^ may, 
however, be a complex number). After simpUfying the curl 
equations in we have the following relationships between 
the different field components 



\u)[i{x)}iy{x) =ikzE^{x) + —E^{x) 



dx- 

ik^Hy{x) =iwe{x)Ex{x) 
_d 
dx 



(4) 



-Hy{x) =io)e{x)Ez{x). 



Using these equations we get the following differential equa- 
tion for Hy 



e{x) 



1 



dx e{x) dx 



x)e{x)^Hy = klHy 



(5) 



and since Ez{0) = by the pec wall at x = 0, the boundary 
condition for Hy under Q becomes -j^Hy{x)\^_^ = 0. The 
equation Q is in the Sturm-Liouville form (Ref jfij ch. 5). 



A. Point Spectrum 

The standard approach (Ref ^ Chap. 5) in the calcula- 
tion of the point spectrum of a Sturm-Liouville equation as in 
^ starts with a redefinition of the space to iL^ — the set of all 
weighted Lebesgue square-integrable functions such that 



e{x) 



dx < 



which implies that the boundary condition at infinity should 
be lim^t^oo Hy{x) = 0. The inner-product in L^, denoted by 
(( •!•))£, is then defined as 



{{f\g)h= j r{x)g{x)^^dx 



(6) 



In order to have a definite metric for L^, the inner product 
should be such that ((/|/))t- > for all / + so that the norm 
of any non-zero vector will be a positive quantity. This in turn 
implies that to have a definite metric, e{x) should be a real and 
positive number for all x. Within the Hilbert space obtained 
by our choice of the inner product ((•!•)) 5, we can write ([5} as 
LHy = k^Hy. The operator 



L = e{x) 



1 



dx e{x) dx 



+ 0)^ fi(x)e{x) 



is self-adjoint since {{f\\-g))e = {{^f\g))e for all/ and g as long 
as e{x) 6 M and e{x) > 0. One can then easily prove that 
the point spectrum of L is purely real (Ref [35} p. 50). The 
lossless dielectric slab waveguide, which satisfies all the criteria 
we mentioned, therefore has a purely real point spectrum. 

Unfortunately, the arguments above fail for the mim wave- 
guide system since the condition e{x) > is no longer 
satisfied^'. The dielectric constants of metals can have negative 
real parts at some frequencies (e.g., in the infrared and visible 
regions), and generally also have imaginary components cor- 
responding to loss, especially at optical frequencies. We will 
now separately analyze the lossy and the lossless metal cases. 



1. Lossless Case 

Since |Im(em)| « |Re(em)| at infrared frequencies, it is 
worthwhile investigating the case of real, negative permittiv- 
ity, i.e., = -|e„, |. The standard Sturm-Liouville theory is 
not appUcable in this case, because it requires the weighting 
function e{x) to have the same sign over its entire domain of 
definition (Ref [35] p. 50). However, > whereas < 
for the MIM waveguide, under the approximation of negligible 
loss. The definition of the inner- product ([6} becomes indefinite 
in this case, since we can have < for some / ^ 0. As a 

result, we no longer can operate in the Hilbert space. The space 
of functions with an indefinite metric is called the Krein space. 
In contrast to the Hilbert space case, the spectrum of the self- 
adjoint operators in Krein spaces is, in general, not real (Ref 
I36! p. 220). An early analysis of a real Sturm-Liouville equation 
with a complex point spectrum can be found in Ref |37] 

To prove that Q accepts complex solutions even when 
e{x) e M, let us work in the well defined space with an 
inner-product as defined in Because (-I-) is always definite, 
we are back in the Hilbert space, but L is no longer self-adjoint 
in Two integrations by parts'* give as 

A 1 A 

e{x) + u)^ fi{x)e{x) 



dx e{x) dx 
with boundary conditions 

Hy{x)\ =0 

^£(x)Hy{x)\ =0. 
dx ix=o 

We see that = e"' Le which makes L by definition pseudo- 
Hermitiai^SE2\ It has been proved that a pseudo-Hermitian 
operator does not have a real spectrum if £^ is indefinite (Ref 
^ Th. 3). 

Alternatively, we can approach the problem by defining 

1^, _ ^ ^ ^2 ^^-j 

dx e(x) dx 
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and rewriting Q as 



fc2 

L Hy = Hy. 

e 



Using iQ it can be shown that L' is self-adjoint in so that 

(H^2|L'H^i) = {\-'Hy2\Hy,) 

which leads to the generalized eigenvalue problem for the self- 
adjoint operators L' and e"^ as 



L'Hv 



where is the eigenvalue. The point spectrum of the self- 
adjoint generalized eigenvalue problem will be complex only 
if both L' and e^^ are indefinite (Ref "41? p. 38). The indef- 
initeness of e"^ is trivial because epsilon can be a positive or 
negative quantity now. To show that L' is indefinite, observe 
that by using the boundary conditions in integration by parts, 
one can get 



e{x) 



dx 



Ax 



which can be positive or negative depending on the choice of 
Hy{x) € L^. Therefore, L' is indefinite, and kl will accept com- 
plex values. Note that the classification of the point spectrum 
into the real and the complex categories is based on kl and not 
kz- Hence, the set of modes with purely real and negative fc^ — 
which leads to a purely imaginary k^ — are categorized as real 
modes in this approach. 



2. Lossy Case 

As we mentioned earlier, £„ has an imaginary part. As a 
result, for those cases when neglecting the imaginary part of 

is not desired, L cannot be made self-adjoint by a redefini- 
tion of the inner-product. Therefore, the point spectrum — the 
set of kl for which (L - k^) does not have an inverse — will 
be complex. A general classification of the spectrum for non- 
self-adjoint operators is still an open problem (Ref p. 301 
and Ref 42). Also, completeness of the spectrum is difficult to 
prove. However, the mim waveguiding problem can be shown 
to have a spectrum which forms a complete basis set even when 
L is non-self-adjoint (Ref ,31, pp. 333-334, Th. 5.3). 



3. Mode Shape 

The dispersion equation that should be solved in order 
to find the k^ values for the modes of the mim waveguide 
is derived by satisfying the continuity of tangential electric 
and magnetic fields at material boundaries and applying the 
boundary conditions. We refer the reader to Ref 31, pp. 462- 
470 and Refs. 8 11 14 18 22 23 for the details. The eigenvectors 
(fn) and the dispersion equation for the corresponding eigen- 



values {kl „) of 1^ for the even tm modes of the mim wave- 
guide are 



(cosh(K,_„A:) 
cosh(Ki,„a) 
-K,„,„(X-£I) 



< X < a 
a < X < 00 



tanh(K,,„fl) 



(7) 

(8) 
(9) 



where Re(Km „) > so that f„{x) does not diverge and is in- 
tegrable. Here « is a discrete index for the eigenvalues and the 
eigenfunctions. Note that we have chosen to write the modal 
shape in terms of the surface mode formulation of Ref [ul 
Sec. III. Surface modes are the main propagating modes for 
the MIM waveguide and have hyperbolic modal shapes. It is 
equivalently possible to describe the modes in terms of oscilla- 
tory shapes using trigonometric functions" — analogous to the 
modes of the dielectric waveguide. Analytical continuation of 
the modal parameters (k, , , /c^ ) makes the two formulations 
equivalent. 



B. Continuous Spectrum 

In this section, we will mathematically show how a contin- 
uous spectrum can exist^^ in the mim waveguide and relate 
it to the continuous spectrum of the dielectric slab waveguide. 
The utility of the continuous spectrum will be evident in the 
mode matching analysis. 

As clearly argued in Ref [43) p. 16, the condition of square 
integrability of the modes can be replaced by the weaker con- 
dition of finiteness of the modes in their domain of defini- 
tion. For the mim waveguide, this would imply a non-zero, 
yet finite electromagnetic field at infinity. These infinite-extent 
and, therefore, infinite energy, continuum modes (which can 
be normalized through the use of the Dirac delta distributions 
as illustrated in Ref 44, pp. 134-135 and Ref 45 pp. 141-148) 
are integrated to realize any physically possible finite energy 
field configuration. In this respect, such an approach is simi- 
lar to the well-known Fourier transform methods, where finite 
energy functions are expanded in terms of the infinite energy 
exponentials. 

Constraining fields to be finite, instead of zero, at infinity 
leads to the following field profile 



(j()v(x) = Ho 
C 



' COsh{Ki^yX) 

cosh(K,_vfl) 

cosh (Km, - a)) 

-I- Csinh(K:m,v(^ - «)) 



< X < a 



a < X < oQ 



■ tanh(K;i,vfl) 



(10) 

(11) 
(12) 



which is calculated very similarly to the dielectric slab example 
in Ref |3ij pp. 462-470. Here v is a continuous index for differ- 
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FIG. 4: (Color online) Spectrum of the mim waveguide for e„ = -143.497 and 2a = A/4 where A = 1550nm is the wavelength of operation. 
There are four real modes and an infinite number of complex modes, all denoted with the • symbol. The thick line denotes the continuous 
spectrum. Due to the fact that e,„ is real, complex modes come in complex conjugate pairs. Insets show the Hy mode shapes in the x direction 
for the discrete spectrum (tmo through TMg) and the continuous spectrum (tmci and tmcz) — solid lines in the insets are the real part of the 
mode, dashed lines are the imaginary part. The locations of the drawn continuous modes are shown by the • symbol. Modes in the continuous 
spectrum are purely oscillatory in the x direction. Complex modes have a small decay, which is not visually apparent in the inset for tms . Arrows 
(t) denote the position of the modes of a parallel plate waveguide with a separation of a — equivalent to the limiting case h ^ O'm Fig.ls] 



ent functions in the continuous spectrum. For finite 0v> the ar- 
guments inside the hyperbolic functions for x > a in ([loj, Km,v, 
should be purely imaginary which implies that Re{ic^^^) < 
and Im(K^^) = 0. These conditions can be written in terms of 
kz^v by using as 

Re{kl^^ - (1)^ l^em) < and Im(A:^ „ - w^^em) = 0. 

Note that when dsj holds true, we have ( = -1 in which 
makes ([loj and Q equivalent. 



C. Residual Spectrum 

We saw that is a second order differential equation which 
could also be written as LHy = k^Hy where L is a differential 
operator. In Ref [33] p. 200 it is claimed that differential op- 
erators have an empty residual spectrum, but a proof is not 
provided. In Ref 31 p. 224 the residual spectrum is said not to 



occur in most of the electromagnetic applications, and in Ref 
31 p. 238 the residual spectrum is said to be empty for typi- 
cal differential operators, though it is highlighted that such a 
fact is not a general result. For Q we did not find any vectors 
belonging to the residual spectrum of L. 



D. Orthogonality relationships 

Orthogonality and completeness are two very valuable prop- 
erties of modes, which make the mode matching technique, to 
be discussed in the next section, possible. OrthogonaUty of the 
modes for the self-adjoint, and non-self-adjoint cases are usu- 
ally expressed using different definitions of the inner product. 
In this work, we will use the pseudo-inner product, [•[], defined 
as (Ref 41, p. 42) 

f{x)g{x)dx. 



MK/K) 




FIG. 5: (Color online) Spectrum of the mim waveguide as a function of d = 2a, Re(fcz/fco) and Im(fcz/fco) where ko = 2n/X is the free space 
wave vector. e,„ = -143.497 is assumed, which corresponds to the real part of the permittivity of silver at a wavelength of A = 1550nm. The thick 
line denotes the Umit of the real spectrum. denote the point spectrum for the d = 2a = \/4 case as also shown in Fig. [4] It can be seen 

that only the lowest order mode, tmo, is propagating for d = A/4 and the rest of the modes are highly evanescent. 



It can be shown that two different eigenfunctions of L, 
and fiix), corresponding to two different eigenvalues ^ and 
fc^ 2 are pseudo-orthogonal with e^^{x) weight (Ref [30lP- 330 
and Ref^p. 47) 

[e^VilVi] = 0. (13) 

From Q it can be seen that Vi proportional to the trans- 
verse electric field component of the mode. Therefore, the 



orthogonality condition can also be written as 

E,i{x)Hy2{x)dx= y^Ei(r) xH2(r)- dA. = 

which is the well known modal orthogonality condition proved 
by the Lorentz reciprocity theorem (Ref 46 p. 336), where A 
denotes the cross section of the waveguide. 

One can directly verify ([13} by integration and using ^ - 
2 ~ ^1 2 which is a result of The following orthog- 
onality conditions between the elements of the point {f„) and 
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the continuous spectrum can similarly be proved 

[e"V«|'i''v] =0 for all n and v, 
[e"Vf<|<i''v] =0 for v + ^. 

The orthogonality conditions talked about in this section 
can also be described in terms of the bi-orthogonal relation- 
ships between the eigenfunctions of the operators L and as 
has been done in Ref |3i| Sec. 5.3 and Ref [47] In Ref [48] four 
examples which illustrate how to choose the weight of the in- 
ner product definition so as to have orthogonal basis functions 
are given. 

In the following sections, we will be working with fields at 
the junction of two different waveguides. For notational abbre- 
viation we will use the following convention 



"{i.j;} 



j{I,J!} 



H 



{L.R} 



where {L, R} is used to denote the modes of the left and right 
side of the junction, which leads to the following orthogonality 
condition 



['^{L,R}\'^{L.R} 



L,R} 



(14) 



where 5,j is the Kronecker delta function and fl is the overlap 
integral of the electric and magnetic transverse fields. 

After the classification and analysis of the mim waveguide 
modes, we will now visuahze different parts of its spectrum 
by finding the zeros of the respective dispersion equations 
through the use of the argument principle method as explained 
in Appendix A. We will use the adjectives in Table[l]to further 
differentiate between the modes. 



TABLE I: 
Signifier 
Leaky 
Proper 
Improper 
Forward 
Backward 



Adjectives 

Signified 
Re(Km) < 
Re(Km) > 
Re(K„) = 
Re(fcz) > 
Re(fcz) < 



Leaky modes are not normaUzable and are not part of the 
spectrum. Proper modes can be normalized by the usual in- 
tegration and they form the point spectrum. Improper modes 
can be normalized by using the Dirac delta functions, S{x). 
They form the continuous spectrum. Forward modes have a 
positive phase velocity, whereas the backward modes have a 
negative phase velocity. We decide on the sign of Re(A:2) based 
on Im(fcz): By definition, all modes are propagating in the +z 
direction. Therefore, in the limit z ^ 00, the fields should go 
to zero. Such a behavior is possible only if Im(fcz) is negative, 
since the fields have an exp{-ikzz) dependence. The argument 
principle method gives us the Km value for the modes. By using 
([9} we get the k^ value. We then calculate (k^y/^ and choose 
the root which satisfies Im(fc2) < 0. Different definitions of 



forward and backward modes — including the ones we use — 
are analyzed in Ref ,49, The definition we use for the leaky 
modes is the same as the one used in Ref [5} Sec. 1.5. 

In Fig. Hthe spectrum of an idealized lossless silver-like mim 
waveguide is shown on the plane of for e„ = -143.497 
which is the real part of the permittivity of silver at a wave- 
length A of i55onm (Ref |50|5i| i. There are four real modes 
for 2a = A/4 — tmq, TM2, TM4, TMe — indexed according to the 
number of zero crossings in Hy. There is also an infinite num- 
ber of complex modes, which are those with eight and more 
zero crossings in the insulator region. These modes have a 
with a positive real part that is rather small compared to the 
imaginary part — this can also be deduced from the scale of the 
imaginary axis of Fig. [4 The continuous spectrum is illus- 
trated by the thick line which corresponds to Re(K^) < and 
Im(K;^ ) = 0. This line is also the branch cut of the square root 
function that is used to get from k^. 

The field profiles of the modes in the insulator region, as 
shown in the insets of Fig. [4j look quite similar to the field pro- 
file of the even modes of a parallel plate waveguide with a plate 
separation of 2a. The even modes of a parallel plate waveguide 
have „ = -n^n^ I a^ . We plotted the corresponding „ val- 
ues on Fig. [4]by using ([9} . It can be seen that such a description 
gives a quite good estimate of the location of the modes on the 
complex plane. 

K^, = is the bifurcation point for the point spectrum when 
is purely real. For positive k^, the point spectrum has 
real modes, whereas for negative k^, the point spectrum spHts 
into two branches that are complex conjugates of one another. 

= corresponds to k\ = Cmkl which then impHes k^ = 
-i^\em\ko — bounded modes should have Im(fcz) < 0. In Fig. 
[sjthe point spectrum is visualized as a function of the insulator 
thickness d = 2a of the mim waveguide and the real and imag- 
inary parts of k^ for the same as in Fig. [4] The thick line in 
Fig. [sjdenotes the limit of the real point spectrum which is the 
bifurcation point of the modes. For Im(fcz) < -^/je^fco there 
are no real modes. The point spectrum for the d = 2a = A/4 
case of Fig. ^is also highlighted with spheres in Fig. [5] For the 
sake of clarity, we only drew one branch of modes after the bi- 
furcation. Due to the scale of the axes in Fig. [5] the quite small 
real part of k^ for the complex modes after the bifurcation line 
is not visibly discernible, but is numerically there. 

In Fig. [6]the spectrum is plotted this time on the plane 
with the same set of parameters as used in Fig. [4] The two 
branches after the bifurcation form the forward and backward 
proper, complex modes. Leaky modes which are not part of 
the spectrum, but nevertheless are solutions to l|8]l-l|9} are also 
shown. When loss is introduced to the metal, the spectrum 
moves on the complex plane as illustrated in Fig. [7] Forward, 
proper, complex modes of Fig.[6]turn into leaky modes by mi- 
grating into the third quadrant of the complex plane. 

In Tablejnjwe provide the numerical values of K„/ko, where 
fco = 2n/X, for the modes in the point spectrum as labeled in 
Fig. [4] The upper line in each row is the value for the e K 
case, the lower line is for the e C case. The effect of loss is 
greatest on the TMe mode. 
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FIG. 6: (Color online) Visualization of the spectrum on the k„ plane. 
All parameters are the same as in Fig. [4] Proper and modes are 
shown for forward, •, and backward, ©, modes. The thick line is the 
continuous spectrum. The complex modes have a small, yet non-zero, 
positive real part — which is what makes them proper modes, unlike 
the continuous part of the spectrum which is improper Thus complex 
modes and the continuous spectrum do not intersect. 



TABLE II: K,„/ko values for the modes labeled in Fig. [4] 



TMo 


12.02521374394057 

12.03170325421919 + i 0.39535643587620 


TM2 


11.34454132059978 

11.35226216175467 + i 0.41892103144838 


TM4 


8.98087712606770 

8.99644754875734 + i 0.52888321538511 


TMs 


0.7136870643968289 

2.247924588647662 + i 2.124681976891650 


TMg 


0.00887301858491 + i 10.55951095636977 
0.45907739584359 + i 10.56851955358188 



IV. MODE-MATCHING 



In this section, we will make use of the spectrum of the 
MiM waveguide to calculate the scattering at the junction of 
two guides with different cross sections. We will use the mode 
matching techniqu e comm only used in the microwave and 
the optical domain&sMUllIU 



A. Birth of the Discretuum 

The presence of a continuous spectrum leads to the forma- 
tion of integral equations when the mode-matching method 



FIG. 7: (Color online) Effects of loss on the spectrum. All parameters 
are the same as in Fig. [ejexcept e„, = -143.497 - 19.517. Comparison 
with Fig. [e] shows that the forward, proper, complex modes on the 
fourth quadrant have moved to the third quadrant and thus became 
leaky modes. 



is applied (Ref. [44) Ch. 5). The integral equation is then ex- 
panded using an orthogonal basis set — not necessarily that of 
the modes — to solve the scattering problem. 

Another way to approach the scattering problem is to limit 
the transverse coordinates by a pec wall. This approach has 
the effect of discretizing the continuum part of the spectrum 
(Ref. [43] pp. 38-41 and Ref. I36] pp. 204-205) — turning it 
into a discretuum^^ . To limit parasitic reflections from the pec 
walls, absorbing layers can be positioned before the pec termi- 
nation (Ref. 53, Ch. 3). In Ref. [59] Sec. 3.2b detailed analysis of 
how the continuous spectrum appears from a discrete collec- 
tion can be found. We will use a pec wall to discretize the con- 
tinuous spectrum. Also, we will not use any perfectly matched 
layers to limit parasitic reflections since the metallic sections 
with permittivity effectively absorb the fields away from the 
junction. 

The geometry is as shown in Fig. [s] For the left waveguide 
the dispersion equation for modes becomes 



tanh(K:,_„a) = m.w/Sm ^anh(Km,Hfe) 



(15) 



which asymptotes to ([sj as h 00. The transverse magnetic 
field shape is 



f„(x) = Ho 



'cosh(K,_„x) 

cosh(K,_„fl) 

cosh {Kfn,n(x - a- h)) 

COSh{Km,nh) 



Q <x < a 



a < X < a + h 



(16) 
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In Fig. [9 we show the effects of positioning a pec wall at 
the top of the lossy mim waveguide. The point spectrum is 
almost the same as in the case without a pec wall at the top. 
The continuous spectrum is discretized, and shows an anti- 
crossing behavior (the repulsion between modes which couple 
to each other) similar to the one observed in coupled wave- 
guide systems. One way to understand the anti-crossing is to 
get rid of the pec walls by the method of images to come up 
with an infinite lattice of parallel mim waveguides each sepa- 
rated from each other by a distance 2{a + h). We observed that 
the perturbations to the discretuum decrease as we increase h, 
as expected from coupled mode theory. Also, the magnitude of 
anti-crossing behavior in the discretuum depends on the dis- 
tance to the nearest mode in the point spectrum. The closer 
the point spectrum gets to the continuous spectrum, the larger 
the anti-crossing effect is. We should note that the modes of 
the PEC terminated mim waveguide are equivalent to the even 
modes of a one-dimensional metalUc photonic crystal at the 
center of the Brillouin zone as described in Ref 23 Sec. iv. 
More information on the evolution of the modes and their de- 
pendence on material properties can be found there. 




1/200 2/200 3/200 
Xw) ' ' 



® 



1 

Re(Kj/fc„ 



B. Are the modes complete? 

Before we attempt the calculation of the scattering proper- 
ties of modes at waveguide junctions, we will first investigate 
the completeness properties of the set of modes we have at our 
disposal — the point and the continuous parts of the spectrum. 
The way we will test completeness is to expand the main mode 
of an MIM waveguide of a given thickness in terms of the modes 
of the MIM waveguide with a different thickness. 

Suppose that we are working with the geometry depicted in 
Fig. M Let us expand the fc**^ mode on the right hand side of 
the junction in terms of the modes on the left as 



FIG. 9: (Color online) Effects of termination of the metal by a pec 
wall — as in Fig.jsj— on the spectrum. e,„ = -143.497 - i9.517, e, = 1.0, 
2a = A/4, h = 2A, ko = 2n/X, X = 1550nm. The spectrum looks very 
similar to Fig. [7] In order to highlight the discretization of the con- 
tinuum, we show a zoomed portion of the k,„ plane. The complex 
modes and the discretuum show anti-crossing behavior. We use two 
different axes for the real part of k,„- The scale on the bottom refers to 
the complex modes, ©, whereas the one on the top refers to the dis- 
cretuum, •. The inset shows a portion of the infinite lattice created by 
the repeated application of image reflection off of the pec boundaries. 
Dashed lines in the inset signify the locations of zero tangential elec- 
tric field, where a pec termination can be applied. In separate calcu- 
lations not shown in this figure, we observe that as h is increased, the 
magnitude of the anti-crossing behavior of the continuum decreases 
as a result of decreased coupling between adjacent mim waveguides. 



m=l 

In order to find A(;„ we test the above equation (i.e. discretize 
the equation by the use of integration of both sides by a given 



function) with h^"^ and use ([14} to get 



n 



in) 




FIG. 8: (Color online) Geometry for mode matching, x = a + h plane 
of Fig. [sjis terminated by a perfect electric conductor which leads to 
a discretization of the continuous spectrum. 



Very similarly, we get the following for the magnetic fields 



1 a 



(m) 



(17) 



What is the error in this expansion? We can get a measure of 
it by writing it as 



1 n 



(m) 



Calculating the pseudo-inner product of the above expression 
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with /i^*^' and then dividing it by O^*^' gives an error estimate 



as 



1-E 



Calculating the error based on the magnetic field expansion 
results in the same expression. 

In Refs. I60I61I62I63I64I65I the importance of the complex 
modes has been demonstrated. In Fig. [lojwe show the impor- 
tance of the continuous spectrum. As shown in Fig. [io|a) with- 
out the continuous spectrum, the field expansion converges, 
but to a field profile which is not the same as the desired pro- 
file of the right junction. On the other hand, inclusion of the 
continuous spectrum through the discretization of the contin- 
uum by a PEC wall leads to the correct field profile as illustrated 
in Fig. [io|b). It can be observed that the expansion based on 
the point spectrum only quite nicely fits the field profile of the 
right waveguide in the insulating region of the left waveguide 
{xjX < 0.125); however, in the metal region {xjX > 0.125) the 
expansion fails. The point spectrum of the left waveguide has 
an exponentially decaying field profile for xjX > 0.125 which 
turns out to be insufficient for the expansion of an arbitrary 
field profile in the metal. The continuous spectrum, with its 
non-decaying field profile, makes field expansion in the metal 
region possible. 




C. Field Stitching 

Now that we know how to treat the continuous spectrum 
and are confident that the collection of the point and the con- 
tinuous spectrum results in a complete basis set, we can pro- 
ceed with the mode-matching formalism. We will begin by as- 
suming that the p'*^ mode of the left waveguide propagates to- 
ward the right, scatters and creates the following set of fields at 
the right and left sides of the junction, which by the continuity 
of the tangential magnetic and electric fields, are set equal 



k=\ 



Y, {Smp - Rmp) &i {x) = Tkpe\ (x). (19) 

m=l k=l 

Here Rmp is the reflection coefficient of the m''^ mode of the left 
waveguide in response to an incoming field in the mode. 
Similarly, Ti^p is the transmission coefficient of the k'^ mode 
of the right waveguide. Note that we chose Rmp to denote the 
reflection coefficient for the transverse magnetic fields, which 
automatically results in -Rmp as the reflection coefficient for 
the transverse electric fields. 

In Ref [66l it is shown that the testing of the above equations 
should be done by the magnetic field of the larger waveguide 
for enforcing electric field continuity ([igj and by the electric 
field of the smaller waveguide to enforce the magnetic field 
continuity l|i8]l. Although that analysis was specifically done 



FIG. 10: (Color online) Effects of the continuous spectrum on mode 
expansion. £„ = -143.497 - i9.517. The main mode of the right wave- 
guide, ftjj'^ (x), with 2a' = A/20 is expanded in terms of the modes of 
the left waveguide, /t['"^(x), with 2a = A/4. The dashed lines show 
the magnitude of the magnetic field profile of the main mode of the 
right waveguide. The solid line shows the result of the expansion of 
( [17) . Shaded regions show the error from the use of a particular ex- 
pansion, i.e., the difference between the expansion and the actual right 
waveguide mode. Due to symmetry, only half of the fuU modal pro- 
file is plotted, (a) Open waveguide of Fig. 3] expansion made using 
the point spectrum only, (b) pec terminated mim waveguide of Fig. [s] 
with a + h = a' + h' = X, expansion made using the discretuum and 
the point spectrum. 



for waveguides with perfect metals {\em\ 00), we still use 
that strategy so that the formulation hmits to the correct one 
should the metals be made perfect. 

For those cases where a < a', we will take the pseudo-inner 

product of ([18} with e["^ and of (|i9j with ^jjj"^ Furthermore, 
assuming there are L modes on the left and _R modes on the 
right, we get 

L R 
E (^Smp + Rmp^ Smn — ^ Ti^p 6^ \hj^ 

m=l k=l 
L R 

(^mp ~ Rmp) ^ l^ji ^ j ~ '^kp^^i ^ ^kn 



k=l 
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with the help of ([14}. When a > a' hj using 4"' and h["^ to 
test l[i8]l and ([19} respectively, we arrive at the following set of 
equations 



m=l L I J j^^j 



„p and Tip as the un- 



These are linear matrix equations with R„ 
knowns. After calculating the inner products, the set of equa- 
tions can be inverted to give the reflection and transmission 
coefficients for the modes. 

In Fig. [uj we compare the mode-matching method with the 
finite-difference frequency-domain (fdfd) technique*^. In Ref 
|68 scattering at mim junctions was investigated using fdfd. 
It takes relatively few modes for the mode-matching calcula- 
tions to converge. Without the continuous spectrum, the mode 
matching results converge to the wrong result. Inclusion of 
the continuous spectrum decreases the error to around 2%, 
which is probably due to the space discretization of fdfd sim- 
ulations as well as the method used in the de-embedding of 
the scattering coefficients from fields. As is also evident from 
Fig. ujthe utility of the single mode {L = R = I) mode- 
matching calculations increases as the dimensions of the wave- 
guides decrease. The single mode approximation is closely re- 
lated to the simpUfied impedance model investigated in Ref 
|68]where it was shown that for deep subwavelength structures 
impedance models are a good approximation. In Fig. [ujwe also 
show the effect of neglecting the backward modes in the mode- 
matching calculations for the 2a = O.lA case. Backward modes 
are important in this sub-wavelength geometry; however, for 
the wider geometries of the 2a = 0.5A and 2a = 0.9A cases 
we did not observe any increase in the error when backward 
modes were neglected in the mode-matching calculations. 

Analysis of the convergence of the field expansions on both 
sides of a junction is an important criter ion for as sessing the 
validity of the mode -matching techniqud^3z£lZ!IZ3 In Fig. [12] 
we show the magnetic field profile at the junction of two mim 
waveguides. As is evident from the figure, convergence of the 
fields on both sides of the junction is obtained only when the 
continuous spectrum is also taken into consideration. Oth- 
erwise, the fields just on the left and just on the right of the 
junction do not agree with one another, showing one or both 
calculations to be in error. The clear conclusion from this nu- 
merical illustration is that the point spectrum on its own is not 
sufficient to describe the behavior of the waveguide junctions. 
Inclusion of the continuous spectrum is essential. 

In Fig. [isjwe visualize the scattering coefficient of the main 
mode of the mim waveguide. We do the calculations in two 
different ways, one using fdfd, and the other using the mode- 
matching technique with the point and the continuous spec- 
trum. When applying mode-matching, we use the a > a' for- 
mulation for ijfi calculations and the a < a' one for R^^^. There 
is a very good match between the results of the two techniques, 
verifying the applicabiUty of the mode-matching method. 
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FIG. 11: (Color online) Convergence study of the reflection coeffi- 
cient, R\i, of the main mode of the left waveguide traveling toward the 
right waveguide for 2a = 0.9A, 2a = 0.5A and 2a = O.lA. a ja = 0.4 
and e,„ = -143.497 - i9.5I7 for aU cases. Dashed lines are for cal- 
culations including the point spectrum only. Solid lines are those 
with both the point and the continuous spectrum. Empty circles de- 
note the calculations done with the forward point spectrum and the 
continuous spectrum for the 2a = O.lA case only. Error is defined as 
I (R^^"" - R^-'""") /Rl;""° I where mm stands for mode-matching and 
FDFD for finite-difference frequency-domain calculations. The inset 
shows the junction geometry. 



V. DISCUSSION 

We will begin this section by a comparative study of the 
modal structures for the parallel plate, the dielectric slab and 
the MIM waveguides. Our aim will be to frame the mim case 
as a bridge between the dielectric slab and the parallel plate 
waveguides. 

The two dimensional symmetric dielectric slab waveguide — 
where the cladding and the core are composed of lossless, pos- 
itive permittivities — has no proper complex modes (Ref [46j p. 
718) as expected from the regular Sturm-Liouville theory. The 
dispersion equation for the dielectric slab waveguide — which 
is the same as ([s} — does have complex roots. However, all of 
those complex roots correspond to leaky modes that are not a 
part of the spectrum. 

The point spectrum of the dielectric slab consists of 
real propagating modes with pure exponential decay in the 
cladding region and with positive power flux in the direction 
of propagation, z. The number of modes in the point spec- 
trum is a finite quantity. The cutoff condition for the point 
spectrum is obtained by the equality = 0. The continu- 
ous spectrum of the dielectric slab starts just below the cutoff 
condition for the point spectrum where < 0. In this range, 
the modes extend to infinity in the transverse direction, but 
they remain bounded. We can divide the continuous spectrum 
into two sections based on the sign^' of 



+ w^iAo \em I- The section where A:^ > is called the radiative 
part whereas the section fc^ < is called the reactive part of 
the continuous spectrum (Ref [44jPP- 128-132 and Ref 5 pp. 
19-28). For both the radiative and the reactive parts of the con- 
tinuous spectrum, the transverse field profiles of the modes are 
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FIG. 12: (Color online) Magnetic field Hy at the junction of two mim 
waveguides as a result of the scattering of the main mode of the left 
waveguide traveling toward the right — the inset shows the schematic 
of the geometry. The left waveguide has 2a = 0.9A, the right wave- 
guide has 2a' = 0.36A so that a' /a = 0.4 as in Fig. '^a + h = a' + h' = 
3 A/4. Fields on the left of the junction are dashed, fields on the right 
are shown with solid lines for both the real and the imaginary part 
of the magnetic field profile. The diflierence between the left and right 
fields is shaded. Due to symmetry, only half of the field profile is plot- 
ted. Vertical dotted lines at x = 0.18A and x = 0.45A denote the end 
of the insulator region for the right and the left side of the junction 
respectively, (a) Mode-matching calculations using the point and the 
continuous spectrum — 6o modes in total — showing good agreement 
between the fields just on the left and just on the right of the junction, 
(b) Mode-matching calculations using the point spectrum only — loo 
modes in total — showing clear disagreement between the calculated 
fields just on the left and just on the right of the junction. 



sinusoidal standing wave patterns. As the sinusoidal variation 
in the transverse direction (k^ = i|K„,|) becomes more rapid, 
the sinusoidal variation in the propagation direction (fc^) de- 
creases. The jump from the radiative to the reactive part oc- 
curs when the sinusoidal variation in z goes to zero: = 0, 
Km/^o = i\/kml- the reactive part, the variation in x is so 
rapid that, modes decay as they propagate in z. The radia- 
tive part is oiten visualized in the mind's eye by plane waves 
that originate in the cladding region far away from the core, 
propagate at an angle towards the core, reflect off of the core 
and propagate away from it. The interference pattern between 
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FIG. 13: (Color online) Reflection coefficient of the main mode at the 
junction between two mim waveguides of different insulator thick- 
nesses, plotted on the complex plane within the unit circle. Both 
waveguides have e„ = -143.497 - 19.517 and e; = 1.0. Filled circles, •, 
are fdfd results, empty circles, o, are mode-matching results. The left 
waveguide thickness is fixed at 2a = 0.9A. The right waveguide thick- 
ness varies from 2fl' = {0.02A, 0.04A, . . . , 0.9A}. The origin is the 
zero reflection point that corresponds to a = a'. As a' decreases pro- 
gressively toward zero, we move progressively along the curves away 
from the origin. The first set of curves, R^, are for the case when 
the mode of the left waveguide, traveling from left to right, is scattered 
by the junction. The second set of curve , Rfi, are for the case when 
the main mode of the right waveguide, traveling from right to left, is 
scattered by the junction. Insets illustrate the respective cases. Details 
about the fdfd calculations can be found in Re£[68l 



the incoming and the outgoing plane waves leads to a stand- 
ing wave pattern in the transverse direction and a propagat- 
ing plane wave in the positive z direction. The reactive part 
is harder to think about in terms of plane wave propagation 
since in this part of the continuous spectrum is purely imag- 
inary and, therefore, the modes are decaying in the positive z 
direction. The plane wave picture is often extended to the re- 
active range by allowing for the possibility for the plane waves 
to come at an 'imaginary' angle of incidence. 

The spectrum of the parallel-plate waveguide with pec 
boundaries is less nuanced than that of the dielectric slab case. 
There is no continuous part of the spectrum. All modes belong 
to the point spectrum and form a discrete basis set. There are 
infinitely many discrete modes supported by the parallel-plate 
waveguide. Finitely many of them are propagating modes and 
carry a positive energy flux in the z direction. The remaining 
ones do not carry any energy and are called evanescent. 

The spectrum of the mim waveguide is a hybrid of the dielec- 
tric slab and the paraUel-plate waveguides' spectra. The mim 
waveguide has both the point and the continuous spectra. Let 
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us start with the continuous spectrum first. 

In Section iii we have shown that the continuous spectrum 
of the lossless mim waveguide exists for all < w^jAe„ which is 
equivalent to < -w^jA\em \ since is a large negative num- 
ber for the MIM case. We see that the continuous spectrum 
of the MIM waveguide is purely reactive and can be thought 
of as composed of plane waves coming at an imaginary angle 
through the metal and reflecting off of the insulator region. 
Equivalently, at a more fundamental level, they are the solu- 
tions of the wave equation with the condition that the fields be 
finite and non-zero at infinity. 

The point spectrum of the mim waveguide has infinitely 
many members, similar to the parallel-plate case. Indeed, one 
can think of the point spectra of the mim and the parallel-plate 
waveguides as analytical continuations of each other as illus- 



trated in Ref '9' for the main tmq mode of the mim and the 
TEM mode of the parallel-plate waveguides. Consequently, the 
main tmq mode of the mim waveguide can be thought of as the 
symmetrical coupling of the two surface plasmon modes at the 
top and the bottom metal-insulator interfaces^"*. 

There are no evanescent discrete modes in the dielectric slab 
waveguide, but the mim waveguide supports them. As illus- 
trated in Fig. [5] there are infinitely many discrete modes for 
any given insulator thickness. Finitely many of those are real 
modes {k^ e M), the remaining ones are complex (fc^ e C). Of 
the real modes, only those with fc^ e M carry any power flux. 
These observations are strictly true only for the case e K. 
When there is loss in the system, all modes do carry small, yet 
finite, amount of power We have verified these claims by cal- 
culating the power flux as 



P=^Re(| E,{x)H;{x)dxy^-\Hl\xRe- 



k^ sinh(KjHd)/Ki, + sm{Ki,d)/Kij 



a)emRe(Km) wCi cosh{Kif.d) + cos{Kijd) 



using Q, and = k:,r + Ik,, for different modes in the point 
spectrum. For waveguide junctions, it has been argued that 
neglecting complex mode pairs on either side of a waveguide 
junction leads to a discontinuity in the reactive energy stored 
at the junctions'''. A similar argument can also be made for 
the highly evanescent continuous modes of the mim geometry. 
Real and complex bound modes are exponentially decaying 
in the metal region. At the junction between two mim wave- 
guides, the smaller waveguide's discrete bound modes cannot 
account for the field leakage into — and therefore reactive en- 
ergy storage in — the metal region due their exponential decay. 
The continuous modes which extend infinitely into the metal 
region make it possible to account for the leakage into metal 
regions. 

We have shown the necessity to take into account the full 
modal structure of the mim waveguide by the calculations we 
presented in Section iv. Here, we should note that, when we 
did mode-matching calculations for the lossless mim geome- 
try, we occasionally observed convergence problems while we 
were trying to reproduce the fdfd results. However, when we 
included loss, all our calculations converged and we did repro- 
duce the lossy fdfd results as we have illustrated in Fig. [is] The 
matrix that one needs to invert to solve the mode-matching 
equations has a higher condition number in the lossless case 
compared to the lossy one. That may explain the difficulties 
we faced. 

In the remainder of this section, we will draw some connec- 
tions between optics and other branches of science with the 
hope of expanding the analogical toolset we use for analysis — 
as exemplified in Ref [75] 

The one dimensional Schrodinger equation and the electro- 
magnetic wave equations in layered media are closely related. 
Both are in the Sturm-Liouville form and one can map the dis- 
persion equation for the te mode of a dielectric slab waveguide 



to the dispersion equation for the modes of a finite potential 
well (Ref I76] p. 11). Furthermore, if one allows for the fi- 
nite potential well to have different effective masses in the well 
and the barrier regions, then a mapping to the tm mode dis- 
persion equation ^ also becomes possible. In Refs. |77j78|79| 
exact closed form analytical solutions for the modes of a sin- 
gle, finite quantum well are developed. These solutions can be 
mapped to the transverse electric modes of the mim or the di- 
electric slab waveguides and perhaps with some labor could be 
expanded to the tm case as well. Effects of discontinuities in 
quantum well potentials are shown to lead to changes in the 
reflection spectrum of the wells in Ref 80 It is intriguing to 
ask whether such studies could be useful in optics for the in- 
vestigation of the effects of material interfaces. Recently, it was 
shown that non-hermitian potentials in the Schrodinger equa- 
tion can have purely real spectra due to the certain symmetries 
of the Hamiltonian of the system''. Pseudo-hermiticity, which 
we have touched upon in Section iii, has been shown to play 
an important role in the interpretation of these systems*^. In 
Ref [83|a parallel plate waveguide with impedance boundary 
conditions was analyzed by the help of the definition of an in- 
ner product which reveals some hidden symmetries'"' '' of the 
system. The operator theoretic findings summarized in Ref .STl 
can have implications for the analysis of waveguides. 

Lastly, in the microwave literature, the unique definition — 
if there is any — of the impedance of an arbitrary waveguide 
mode is an active area of research. The causal waveguide 
impedance definition of Ref 86 seems to formulate a uni- 
fying framework to merge different interpretations together. 
It seems worthwhile to ask what the tmq mode impedance 
of an MIM waveguide would be for a causal definition of the 
impedance given the Kramers-Kronig relationships for wave- 
guide modes as investigated in Ref |^ 
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VI. CONCLUSION 

In this paper we investigated the even tm modes that the 
MiM waveguide supports. We based our analysis in the lan- 
guage of operators and used methods developed for Sturm- 
Liouville systems to expand the results reported in Ref [23] 
The mathematical structure of the odd modes are very simi- 
lar to the even ones and can be derived in a similar manner. 
These findings were in accordance with Ref 26 where it was 
shown that in general, open structures will have complex and 
continuous spectra. 

After the investigation of the modes, we showed their util- 
ity and relevance by the mode-matching method. We investi- 
gated the problem of modal scattering at the symmetric junc- 
tion of two MIM waveguides with different cross sections''^ and 
successfully applied the mode-matching technique to predict 
the modal reflection coefficients calculated by fuU-field simu- 
lations. Lastly, we commented on some of the possible links 
between the quantum mechanics, optics and microwave liter- 
ature and considered possible research directions. 

The knowledge of the set of orthogonal modes which form 
a complete basis for a given geometry leads to a much more 
simphfied algebra and speeds up calculations. The results 
of this paper are valuable for electromagnetic scattering cal- 
culations involving the mim geometry. Our results would 
also help in the analysis of optics experiments involving mim 
waveguides"* Furthermore, the results reported are also 
useful for analyzing plasmonic quantum optics^3^3111Il and 
Casimir effect deviceS^^l^. 

The analysis made in this paper can be generalized 
for other related g eometries involving metals at optical 
frequencie^SZmillliSSI j^y^ j-j^-j^ Qf modes available in 
the MIM geometry suggests that modal investigation of ar- 
bitrary three-dimensional nano-metallic waveguides — which 
are thought to replace the electrical interconnects on future 
computing devices — will require novel means of deducing 
their discrete and continuous spectra. 
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APPENDIX A: ARGUMENT PRINCIPLE METHOD 

In this section we will briefly describe the method we used 
to find the zeros of ([sj and One of the main of problems 
of any root finding algorithm is the starting point in the search 
domain. When the search needs to be done over two dimen- 
sions of the complex plane, brute force approaches have lim- 
ited appUcabiUty. Luckily, there have been many advances in 
the root finding algorithms for waveguides^"' ioz.io3.io4,i°5 h)6^ 



If the function f{z f^is analytic, possesses no poles on and 
within the closed contour C and finally if f{z) does not go to 
zero on C then: 

where Zk denote the zeros of f{z) in C, N is the total number 
of zeros in C and m is an arbitrary nonnegative integer. Specif- 
ically, for m = one gets the total number of zeros within C. 
Knowing the number of zeros in a closed region enables one to 
dismiss regions of the complex search space in which there are 
no zeros. Furthermore, given the number of zeros in a region, 
one can do a subdivision until there is only one zero in the re- 
gion of interest as described in Ref 103 For those regions with 
a single zero, another contour integral as in l |Ai| l with m = 1 will 
give the location of the zero. If better numerical accuracy is de- 
sired, one can do a local search given the approximate position 
obtained by the contour integration. This method of finding 
the zeros of a function by repeated integration on the complex 
plane is called the argument principle method (apm). 

Note that for the apm to work, the function f{z) should be 
analytic in C. That requires no branch points to exist in the 
closed contour. Our implementation takes the variable 
of interest. Equations ([sjl and are rewritten in terms of only 
Km by using Singularities of tanh(K, a) and of tanh(Km?i) 
are removed by multiplying both sides of l|8} by cosh(K, a ) and 
of ([15) by cosh(K,a) cosh{K,„h). At this stage the equations 
look like: 

BmKi Sinh{Kia) + EfKm cosh(K,fl) = (A2) 
CmKi sinh(K/fl) cosh{Kmh) + 

CiKm cosh(Kifl) smh{Kmh) = 

with Ki = y^K^ - fco(em - £() where kg = 2n/X. The only 
branch points'"* are those caused by the square root function 
in the definition of k,-. Note that the function \/z sinh{\/z) 
is single valued everywhere on the complex plane, since re- 
gardless of the choice for the sign of the square root, the re- 
sult does not change. Similarly cosh(\/z) is single valued due 
to the evenness of the hyperboUc cosine function. Therefore, 
l\A2\ and l\A^\ have no branch points and are analytic in the 
whole complex plane. 

In our implementation, we also did the conformal mapping 
Kfn = e^ to map the complex plane into strips, z = log(Km) = 
log|Km| + iarg(Km) = log|Km| + iarg(Km) + i2;im for any integer 
m. The logarithm function is multi-valued, and repeats itself 
in strips that result from the i2nm term. For instance, when 
we are interested in finding proper modes with Re(Km ) > we 
only need to search the strip -n/2 < Im(z) < n/2. 

Lastly, using computer algebra systems that can do symbolic 
mathematical manipulations, one can easily calculate f'{z) 
from the definition of f{z) given by ( |A2| l or ([As} in addition 
to doing numerical integrations on the complex plane. 
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APPENDIX B: SOME DETAILS OF THE MODE MATCHING 
ALGORITHM 

Once the full set of modes are found using apm, implemen- 
tation of the mode-matching algorithm reduces to the calcula- 
tion of the overlap integrals to build the matrix equation that 
should be solved to get the reflection and transmission coeffi- 
cients, Rkp and Tj^p, of ([i8}-(|i9]l. Since we have the analytical 
solutions for the fields as given in and overlap inte- 
grals and l^l^''^}] can be calculated analytically 
in closed form. The expressions are too long to reproduce here, 



but a computer algebra system can do the analytical manipu- 
lations. Once we have closed form results for the overlap in- 
tegrals, formation of the matrix equation is very quick. The 
matrices have relatively small sizes and the solution of the lin- 
ear matrix equation proceeds quickly. Also note that we nu- 
merically check the validity of the mode orthogonality condi- 
tion of before we start the mode-matching calculations so 
as to check the correctness of the apm implementation. We 
used Mathematica to implement apm and the mode-matching 
method. The source code will be made available on the Inter- 
net under the gnu general public license. 
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